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Abstract 


Radar  backscatter  measurements  were  performed  on  artificially  grown  sea  ice  on  an 
outdoor  pond  at  the  US.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
(CRREL)  during  the  1990  winter  season.  Measurements  of  microwave  backscatter 
from  simulated  sea  ice  demonstrate  that  the  dominant  backscatter  mechanism  is  the 
surface  scatter.  Both  the  co- polarized  and  cross-polarized  measurements  compare 
favorably  with  the  predictions  of  surface  scattering  models  at  two  frequencies.  The 
surface  parameters  used  in  the  models  were  measured  directly  from  photographs  of 
cross  sections  of  the  saline  ice.  Therefore,  since  two  frequencies  of  backscatter 
measurements  were  predicted  accurately  by  surface  scatter  models,  the  predominate 
backscatter  must  be  from  the  surface. 

These  measurements  were  carefully  calibrated  using  a  complex  vector  error 
correction  scheme  to  remove  polarization  coupling  and  other  systematic  effects.  The 
dynamic  range  of  polarimetric  measurements  made  by  scatterometers  is  limited  by 
coupling  and  polarization  leakage  in  the  antenna  system  and  RF  section.  These 
effects  must  be  accounted  for  and  corrected  if  one  is  to  obtain  accurate 
characterization  of  the  polarimetric  response  of  distributed  targets. 

This  demonstrates  the  need  for  careful  calibration  in  experiments  of  this  type, 
especially  when  cross-polarized  backscatter  measurements  are  performed.  The 
systematic  effects  on  the  cross-polarized  backscatter  introduce  significant  error  and 
may  lead  to  misinterpretation  of  the  results.  This  is  a  vital  consideration  in  the  study 
of  radar  backscatter  from  sea  ice  using  fully  polarimetric  systems. 
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0.  Sea  Ice  Remote  Sensing 


Ice-covered  oceans  have  been  the  subject  of  extensive  research  due  to  their 
importance  in  global  climate  change  and  their  hindrance  to  shipping  traffic  for 
northern  ports.  The  sea  ice  extent  in  the  Arctic  is  extremely  dynamic.  At  its 
maximum  in  March  the  Arctic  ocean  is  covered  by  about  14.8  million  km^  of  sea  ice 
and  the  minimum  extent  is  about  7.8  million  km2  in  September  [Parkinson,  et  al., 
1987].  The  variation  in  sea  ice  cover  itself  is  roughly  equivalent  to  an  area  the  size  of 
the  48  contiguous  states  of  the  US. 

In  global  climate  models  the  predicted  warming  of  the  globe  due  to  increase  in 
greenhouse  gases  is  most  dramatic  in  the  polar  regions.  For  example,  the  high 
resolution  models  of  the  Canadian  Climate  Centre,  the  Geophysical  Fluids  Dynamics 
Laboratory  at  Princeton,  and  the  United  Kingdom  Meteorological  Office  all  show 
more  dramatic  temperature  changes  in  the  Arctic  than  for  the  rest  of  the  world 
[Houghton,  1990].  While  all  three  models  predict  a  global  temperature  increase  of 
between  3.5°  and  4°  C,  they  predict  an  increase  in  the  summer  temperature  in  the 
Arctic  of  between  6°  and  16°  C.  All  three  models  are  equilibrium  models  based  on 
the  instantaneous  doubling  of  carbon  dioxide  in  the  atmosphere. 

The  albedo  of  the  surface  of  the  earth  plays  an  important  role  in  the  amount  of  total 
energy  absorbed  by  the  earth.  The  albedo  feedback  effect  concerns  the  amount  of 
incoming  solar  radiation  that  is  reflected  back  into  space  by  the  earth's  surface  and 
atmosphere. 


5 


The  albedo  of  the  ice  affects  the  amount  of  short-wave  radiation  absorbed  at  the 
surface.  If  the  extent  of  sea  ice  decreases  and  is  replaced  by  open  ocean,  which  has  a 
lower  albedo,  there  will  be  more  absorption  at  the  earth’s  surface.  This  warming  will 
lead  to  a  further  retreat  of  sea  ice  and  an  increase  in  the  amount  of  radiation  absorbed. 
This  is  effectively  a  positive  feedback  effect  in  a  global  warming  scenario. 

In  addition  to  being  a  possible  indicator  of  climate  change,  sea  ice  cover  also 
modulates  energy  exchange  between  the  warm  ocean  and  the  cold  atmosphere.  Sea 
ice  cover  acts  as  an  insulator  between  them.  Different  sea  ice  types  may  also  affect 
this  interaction.  Areas  of  pancake  ice,  for  example,  may  allow  for  more  energy 
exchange  between  the  air  and  ocean  than  an  area  containing  a  multiyear  floe. 

In  addition  to  the  active  roles  mentioned  above  sea  ice  also  has  been  postulated  as  a 
measure  of  the  changing  climate.  Both  the  extent  of  sea  ice  over  a  period  of  time  and 
its  thickness  may  be  a  signal  of  the  onset  of  global  climate  change,  much  like  the 
amount  of  CO2  in  the  atmosphere  and  surface  temperatures  are  pointed  to  today  as 
indicators  of  global  wanning. 

Because  of  these  effects  it  is  essential  in  studying  our  planet  that  sea  ice  is  monitored 
on  a  global  scale.  The  only  way  to  achieve  this  is  by  looking  at  the  earth’s  surface  by 
remote  sensing  from  satellites.  The  primary  problem  in  the  Arctic  regions  for  other 
remote  sensing  techniques  (visible,  optical,  infrared)  is  that  the  polar  regions  are  dark 
for  half  of  the  year.  In  the  other  half  of  the  year  cloud  cover  and  fog  reduce  the 
effectiveness  of  these  techniques.  The  only  methods  for  consistently  monitoring  the 
Arctic  through  darkness,  fog  and  cloud  are  by  active  and  passive  microwave  remote 
sensing.  Of  these  two  methods,  the  active  microwave  remote  sensors,  usually 
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Synthetic  Aperture  Radar  (SAR),  have  much  better  spatial  resolution.  Therefore  SAR 
is  more  capable  of  monitoring  finer  structure  and  sea  ice  dynamics  than  passive 
microwave  devices. 

Because  of  the  importance  of  sea  ice  in  the  global  climate  and  in  its  impediment  to 
ship  traffic  in  the  Arctic,  and  the  capability  of  active  microwave  sensors  to 
characterize  sea  ice,  the  radar  backscatter  signatures  of  sea  ice  have  been  of  intense 
interest  over  the  last  two  decades. 

This  report  concentrates  on  the  understanding  of  the  backscatter  mechanisms  for 
simulated  sea  ice.  The  emphasis  of  the  simulated  sea  ice  work  is  to  produce  high 
quality  measurements  to  be  used  in  the  validation  of  scattering  models  for  use  in  the 
study  of  sea  ice.  The  measurements  for  this  study  were  taken  at  the  US  Army  Cold 
Regions  Research  and  Engineering  Laboratory  (CRREL)  during  the  Winter  of  1990. 
The  results  of  this  study  show  that  the  backscatter  coefficient  is  dominated  by  surface 
scatter. 
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1.  The  1990  US  Army  Cold  Regions  Research  and 

Engineering  Laboratory  Experiment  (CRRELEX'90) 

1.0  Introduction 

Controlled  radar  backscatter  measurements  were  performed  by  researchers  from  the 
University  of  Kansas  Radar  Systems  and  Remote  Sensing  Laboratory  (RSL)  at  the 
US  Army  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL)  in  early 
1990  to  investigate  the  scattering  characteristics  of  sea  ice.  Measurements  were 
performed  at  an  outdoor  saline  ice  facility  using  two  network  analyzer-based  stepped- 
frequency  radar  systems.  These  systems  were  developed  at  RSL  and  operated  at  C 
band  (5.3  GHz)  and  Ku  band  (13.4  GHz).  These  systems  allowed  for  measurement  of 
backscatter  in  amplitude  and  phase  at  four  linear  polarizations:  VV,  HH,  VH,  and 
HV.  Saline  ice  with  varying  degrees  of  roughness  was  investigated  with  and  without 
snow  cover  at  angles  ranging  from  0°  to  60°.  The  CRREL  experiments  were 
performed  to  study  the  relation  of  physical  properties  of  simulated  sea  ice  to  the 
measured  backscatter  under  somewhat  controlled  laboratory  conditions.  The 
simulated  sea  ice  grown  at  the  CRREL  facility  has  physical  properties  which  closely 
model  thin  Arctic  ice  grown  under  quiescent  conditions,  [Swift,  et  al„  1992].  The 
physical  properties  of  simulated  sea  ice  have  been  shown  to  closely  resemble  sea  ice 
found  in  the  Arctic  [Arcone,  et  al.,  1986].  Since  the  ice  properties  can  be  monitored 
closely  in  this  type  of  situation,  the  physics  of  microwave  backscatter  from  saline  ice 
can  be  studied  in  greater  detail  than  in  field  expeditions  where  there  is  no  control  over 
the  environment. 
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These  measurements,  like  all  radar  measurements,  are  limited  by  the  polarization 
purity  of  the  antenna  and  RF  section  as  well  as  coupling  between  transmit  and  receive 
antennas.  The  accuracy  of  cross-polarized  measurements  is  especially  vulnerable  to 
these  systematic  effects  imposed  by  the  radar  system.  Recendy,  emphasis  has  been 
placed  on  utilizing  the  scattering  information  from  all  four  linear  polarizations  in 
studying  the  nature  of  scattering  mechanisms  [Kwok,  et  al.,  1991],  [Drinkwater,  et  al., 
1991].  In  particular,  the  polarization  ratio,  usually  defined  as  the  ratio  of  co¬ 
polarized  o°  divided  by  cross-polarized  o°,  is  often  used  to  categorize  sea  ice.  Also 
the  polarization  “signature”,  or  more  correctly,  polarization  response,  has  been  used 
recently  in  the  study  of  radar  backscatter  from  ice.  This  requires  information  from  all 
four  linear  polarizations,  thus  the  need  for  accurate  cross-polarized  measurements. 

To  obtain  high  quality  measurements  for  comparison  with  theory,  the  calibration  of 
the  scatterometers  must  take  these  effects  into  account.  To  accomplish  this  a  complex 
vector  error  calibration  was  performed  based  on  techniques  described  by  [Riegger,  et 
al.  1989],  [Wiesbeck  and  Riegger,  1991],  and  [Chen,  et  al.,  1991],  but  applied  to 
returns  from  a  distributed  target  (artificial  sea  ice)  as  opposed  to  simple  targets. 

This  technique  allows  one  to  quantitatively  reduce  the  effects  of  polarization  leakage 
of  the  antennas  to  improve  the  dynamic  range  of  the  scatterometer  data.  This 
technique  was  applied  to  both  C-band  and  Ku-band  data  taken  during  CRRELEX'90. 
The  effectiveness  of  the  technique  is  verified  by  applying  it  to  several  calibration 
targets.  These  results  show  that  the  vector  calibration  leads  to  at  least  a  15  dB 
improvement  in  the  dynamic  range  of  the  backscatter  measurements.  The  vector 
calibration  was  then  applied  to  returns  from  artificial  sea  ice.  The  resulting 
measurements  of  scattering  coefficient  from  saline  ice  show  a  much  greater  difference 
between  co- polarized  and  cross-polarized  returns  than  for  the  usual  scalar  calibration 
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technique  [Ulaby,  et  al.,  1982]  in  which  only  a  single  non-depolarizing  target  such  as 
the  Luneburg  lens  or  metal  sphere  is  used. 

After  careful  calibration  to  remove  system  effects,  the  backscatter  measurements  were 
compared  to  the  small  perturbation  model  for  scattering  from  a  random  rough  surface. 
These  results  demonstrate  that  the  experimental  values  correspond  fairly  closely  with 
the  scattering  theory.  In  particular  the  levels  of  the  cross-polarized  measurements  are 
significandy  lower  than  measurements  of  the  cross-polarized  backscatter  that  are  not 
corrected  for  systematic  effects  of  the  radar  system.  The  ratio  of  like-polarized  return 
to  the  cross-polarized  return  also  increases  when  systematic  errors  are  removed. 

In  this  report  the  results  of  the  scatterometer  measurements  are  presented  for  both  the 
case  before  and  after  the  systematic  calibration  is  applied.  The  details  of  the 
calibration  technique  as  well  as  limitations  due  to  assumptions  made  and  the 
measurement  setup  are  presented.  The  vector  calibrated  results  are  then  compared  to 
rough  surface  scattering  models. 
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2.  Step-Frequency  Radar  Systems  Description 

Two  stepped-ffequency  radar  systems  were  used  to  measure  backscatter  from  an 
outdoor  saline  pond  at  the  CRREL  facility  in  Hanover,  New  Hampshire.  The 
measurements  were  performed  at  C  band  and  Ku  band  with  network-analyzer  based 
systems.  Each  system  contains  an  HP  8753  Network  Analyzer  as  the  core.  The 
intermediate  frequency  (IF)  sections  for  both  systems  are  housed  in  a  single  box  as 
are  the  radio  frequency  (RF)  portions  of  the  system. 

The  RF  sections  provided  the  up-conversion  of  the  network  analyzer  output  from  2 
GHz  center  frequency  (IF)  to  a  center  frequency  of  5.3  GHz  for  the  C-band  system 
and  13.9  GHz  for  the  Ku-band  system.  They  also  performed  the  down-conversion  of 
the  received  signal  to  IF  for  measurement  by  the  network  analyzer.  Figure  2.1  shows 
the  block  diagram  for  both  the  C-band  and  Ku-band  systems,  with  particular  emphasis 
on  the  RF  section.  The  network  analyzer  played  the  role  of  the  signal  generator  and 
processor  for  each  radar  system.  It  was  used  to  coherently  measure  the  amplitude  and 
phase  of  the  return  signal  (i.e.  the  S21  scattering  parameter)  at  IF,  for  a  set  of  discrete 
frequencies.  The  frequency  of  the  network  analyzer  is  stepped  through  1000  MHz  of 
bandwidth  for  the  Ku-band  radar  and  500  MHz  of  bandwidth  for  the  C-band  radar  in 
2.5  MHz  steps.  The  network  analyzer  is  used  to  measure  the  return  independently  at 
each  frequency  and  these  measurements  are  stored  on  disk.  Table  2.1  describes  the 
important  parameters  of  the  systems  used.  Photograph  2.1  (a)  and  (b)  show  the  radar 
system  setup  at  the  lower  pond. 


13 


Figure  2.1.  Block  diagram  describing  both  the  C-band  and 
Ku-band  step-frequency  radar  system. 


The  RF  sections  of  the  C-band  and  Ku-band  radar  systems  are  nearly  identical  and  are 
shown  in  figure  2.1.  The  source  in  the  RF  section  is  a  dielectric  resonator  oscillator 
(DRO).  The  output  of  the  DRO  is  mixed  with  the  IF  source  from  the  network 
analyzer  to  generate  the  RF  radar  signal.  The  frequency  of  the  DRO  is  set  so  that  the 
center  frequency  of  the  RF  signal  is  at  5.3  GHz  for  the  C-band  system  and  13.4  GHz 
for  the  Ku-band  system.  This  is  filtered  and  amplified  before  passing  through  the 
switching  network.  The  switching  network  in  the  RF  section  guides  the  signal  to 
either  the  horizontal  or  vertical  polarization  port  to  the  antennas.  This  allows  for  both 
vertical  and  horizontal  polarizations  to  be  transmitted.  The  backscattered  signal  is 
received  through  the  receive  antenna.  The  switching  network  separates  the  received 
signal  into  either  vertical  or  horizontal  polarization.  The  returned  signal  is  then 
mixed  down  to  IF  for  measurement  with  the  network  analyzer. 
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Table^l^te^^rMuenc^Rada^^sten^arametCTS 


Center  Frequency.  (GHz) 
RF  Bandwidth  (MHz) 
Free-space  Resolution 
(cm) 

Number  of  frequency  steps 
Footprint  in  cm  (2m 
range) 

Transmit  power  (mW) 


C  band 

Ku  band 

5.3 

13.9 

500 

1000 

30 

15 

200 

400 

48x48 

30x30 

10 

10 

The  horizontal  and  vertical  polarization  for  the  C-band  radar  system  are  obtained  by 
using  two  pyramidal  horn  antennas  for  receiving  and  two  for  transmitting.  The 
pyramidal  horn  is  also  called  the  standard  gain  horn  since  it  has  a  wide  beamwidth 
and  its  pattern  is  well-documented  [Balanis,  1982].  These  antennas  provided  a  15° 
two-way  equivalent  beamwidth  as  mounted  during  CRRELEX’90. 

At  Ku  band  the  vertical  and  horizontal  polarizations  are  obtained  with  a  single 
diagonal  horn  antenna  for  transmitting  and  an  identical  antenna  for  receiving.  The 
two  antenna  Ku-band  system  is  obtained  through  the  use  of  ortho-mode  transducers 
for  each  of  these  antennas.  The  diagonal  hom  antenna  relies  on  the  superposition  of 
the  TEoi  and  TEio  modes.  In  the  diagonal  hom  these  modes  generate  a  field  pattern 
that  resembles  the  TEi  i  mode  of  a  circular  waveguide  [Love,  1962].  The  far-field 
patterns  that  result  are  shown  here  in  figures  2.2  and  2.3  (taken  from  Love,  1962). 
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The  principal  plane  patterns  demonstrate  that  the  antenna  has  superb  sidelobe 
performance  as  shown  in  figure  2.2.  Ideally,  there  is  no  cross-polarized  pattern  in  the 
principal  planes.  In  the  45°  plane  the  cross-polarized  pattern  of  figure  2.3  is  seen  to 
be  about  15  dB  below  the  normal  antenna  pattern  near  the  main  lobe  in  the  ideal  case. 
This  is  important  when  considering  cross-polarized  measurements  and  the  effect  of 
finite  polarization  purity  on  these  measurements. 

The  network  analyzer  is  controlled  by  a  microcomputer  via  a  Hewlett  Packard 
Interface  Bus  (HPIB)  expansion  card.  A  computer  program  controls  the  stepping  of 
the  network  analyzer  through  the  frequencies  needed  to  produce  the  bandwidths 
shown  in  table  2.1.  The  program  also  obtains  the  scattering  coefficient  measurement 
from  the  network  analyzer  and  is  used  to  store  the  data  to  disks.  The  switching 
network,  which  is  used  to  switch  between  polarizations,  was  manually  controlled  by 
an  external  control  box.  The  switch  network  was  also  used  to  perform  a  delay-line 
measurement  as  an  internal  calibration.  The  network  analyzer,  switch  control  box, 
and  computer  were  housed  in  a  shed  adjacent  to  the  lower  pond  as  shown  in 
photograph  2.2. 
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Figure  2.2.  Antenna  pattern  for  diagonal  horn  (from  Love,  1962).  This 
shows  the  principal  plane  E-  and  H-plane  normalized  antenna  patterns  for  an 
ideal  diagonal  horn  antenna. 
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Figure  2.3.  Diagonal  horn  antenna  pattern  for  45°  planes  (from  Love.  1962).  This 
shows  the  ideal  45°  plane  patterns  for  co-polarized  and  cross -polarized  normalized 
power.  This  shows  the  level  of  cross-polarized  power  as  being  about  15  dB 
below  the  normal  pattern  near  the  main  beam. 
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Photograph  2.1  (a)  Mounting  of  the  C-hand  and  Ku-band  step-frequency 
radar  systems  a:  CRRKL.  lower  pond  facility  tor  CRRRI.r.X’dO  ( h )  Ciose-up 
vj-w  of  the  R!  ho;-  and  a.ntenn;.  mount  shown:;-  the  tv..  -  manorial  Ku-band 
>-r:  ; •  •  v 1  tie  in;-'  C  banc  nvr.anid;:'  horv  antenna- 
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Photograph  2.2  A  network  analyzer,  eonirol  boxes  and  (  oninau 
were  housed  in  a  shed  adjacent  to  C'RRIiLA  lower  poml  for 
campaign. 
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3.  CRRELEX'90  Experiment  Description 


Radar  backscatter  measurements  from  artificial  sea  ice  were  made  at  CRREL’s  lower 
pond  facility  during  January  and  March,  1990.  The  experiment  utilized  an  outdoor 
frozen  pond  which  was  made  with  saline  water  in  order  to  simulate  sea  ice.  This 
facility  provides  a  test  bed  for  experimentation  under  quasi-laboratory  conditions. 
This  allows  for  more  control  than  expedition- type  endeavors.  The  lower  pond  is 
made  of  concrete  and  is  about  10m  x  10m  x  1.3  m,  with  vertical  sides.  The  facility 
has  a  walkway  on  railroad  tracks  which  served  as  a  mobile  base  on  which  the  antenna 
mount  was  placed.  The  RF  sections  were  also  placed  an  the  antenna  mounts,  while 
the  control  and  data  acquisition  systems  were  housed  in  a  shed  adjacent  to  the  pond. 

Measurements  were  made  at  incident  angles  ranging  from  0°  to  60°,  with  all  four 
linear  polarization  combinations;:  VV,  HH,  VH,  and  HV.  Movement  of  the  antenna 
mount  allowed  for  up  to  four  independent  spots  for  reduction  of  fading  which  occurs 
in  radar  backscatter  measurements.  In  addition,  ice  core  samples  were  collected  by 
CRREL  and  analyzed  to  obtain  the  characteristics  of  the  saline  ice.  Ice 
characterization  data  included  measurements  of  salinity  and  temperature  profiles,  ice- 
structure,  including  micro  structure  from  thin  sections  and  macro  structure  (bubble 
size  and  surface  roughness)  from  thick  sections. 
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3.1  Radar  backscatter  measurements 

Radar  backscatter  measurements  and  ice  property  measurements  were  made  in 
January  and  again  in  March  under  a  variety  of  weather  conditions.  Table  3. 1 
summarizes  the  dates,  times  and  angles  of  which  measurements  were  made,  as  well  as 
the  number  of  independent  spots  used. 


Table  3.1 

CRRELEX'90  Lower  Pond  Data  Summary 

Date 

Time 

Incidence 

Frequency 

Spots 

(EST) 

Angles 

7  Jan 

17:00 

12-55 

C,  K0 

3 

9  Jan 

17:44 

0-55 

C,Ku 

4 

24  Jan 

15:44 

0-40 

Ku 

2 

29  Jan 

10:51 

47-55 

Ku 

2 

5  Mar 

13:27 

3-53 

Ku 

3 

6  Mar 

13:33 

3-53 

C 

3 

8  Mar 

21:03 

0-55 

C,  Ku 

4 

10  Mar 

13:48 

0-55 

C,Ku 

3 

11  Mar 

04:58 

0-55 

C,  Ku 

4 

14  Mar 

04:24 

0-55 

C,  Ku 

4 

28  Mar 

06:37 

0-53 

C,  Ku 

3 

In  addition  to  the  backscatter  measurements  from  saline  ice  the  return  from  several 
calibration  targets  were  measured  with  the  radar  systems.  These  targets  included  a 
Luneburg  lens,  a  metal  sphere,  and  a  dihedral  comer  reflector.  These  measurements 
were  made  so  that  careful  calibration  could  be  used  to  remove  the  systematic  errors 
due  to  polarization  leakage  and  other  effects  which  occur  when  using  real  antennas. 


22 


3.2  In-situ  snow  and  sea  ice  measurements 


During  the  lower  pond  experiment  of  CRRELEX'90  the  air  and  ice  temperatures  were 
taken  periodically,  along  with  ice  surface  salinity  and  snow  depth.  These  are  shown 
in  table  3.2.  From  these  parameters  an  estimate  of  the  surface  dielectric  constant  was 
made  using  past  measurements  from  CRREL  [Arcone,  et  al.,  1986].  From  these  an 
approximate  dielectric  constant  was  selected  at  3.75  -  J0.4.  This  is  based  on  the 
temperature  and  age  of  the  saline  ice  sheets. 


Table  3.2.  In-situ  snow  and  ice  information 

Ice  or  Snow  Property 

Date  /  Time 

Temperature 

-9.5°  C  :  snow  surface 

17:00/7  Jan 

-3.0°  C  :  snow-ice  interface 

17:00/7  Jan 

0°  C  :  Ice  Surface 

13:48  / 10  Mar 

-1.5°  C  :  Ice  Surface 

04:58/11  Mar 

-1.0°  C  :  Ice  Surface 

04:24/ 14  Mar 

Salinity 

12  %o  ±  2  %o  at  snow  /  ice 

interface 

9  Jan. 

58  mm 

7  Jan 

57  mm 

9  Jan 

In  addition,  photographs  were  made  of  sections  cut  from  the  saline  ice  sheet.  Any 
snow,  if  present,  was  removed  from  these  sections.  These  were  used  to  characterize 
the  surface  in  terms  of  rms  roughness  and  correlation  length  for  comparison  to 
theoretical  models.  The  photographs  used  to  find  the  surface  parameters  for  the 
experiment  are  shown  on  the  following  pages.  The  first  six  were  taken  on  20  January 
and  the  last  two  were  taken  in  the  latter  portion  of  the  campaign  on  1 1  March. 
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The  procedure  for  obtaining  the  surface  roughness  was  performed  by  digitizing  the 
ice  section  photographs  by  scanning  into  a  computer.  These  were  then  expanded  by  a 
factor  of  10  using  a  graphics  routine.  From  the  expanded  version,  x-axis  and  y-axis 
values  were  read  from  the  image.  The  spacing  between  samples  was  1  mm  and  the 
precision  of  the  height  reading  was  within  1/10  mm,  on  the  expanded  images.  The 
standard  deviation  of  the  surface  height  (surface  roughness)  was  calculated  from 

(3.1) 

where  N  is  the  number  of  samples  in  the  sequence  and  z  is  the  mean  value  of  the 

surface  height  [Ulaby,  1982].  The  autocorrelation  of  the  surface  height  is  also 

needed  in  scattering  models  and  was  obtained  from, 

n-i 

X(z,  -z)(zM  -z) 

P(X)  =  J-1 -  - ,  (3.2) 

N 

where  N  is  the  number  of  samples,  and  x'  is  the  spatial  displacement  defined  by  the 
sample  number  and  the  spacing  as  x’=(j+l)Ax. 

The  correlation  length  was  obtained  from 

P(LC)  =  -,  (3.3) 

e 

where  Lc  is  the  correlation  length  and  as  usual, 

e  =  lim(l  +  ±)‘ -2.718.  (3.4) 

The  computed  autocorrelations  for  the  slabs  taken  in  the  January  experiment  are 
shown  in  figures  3.1,  3.2,  and  3.3.  Figure  3.1  shows  the  correlation  for  slabs  la  and 
lb,  with  slab  la  having  a  longer  correlation  length.  This  is  due  to  the  step-like  jump 
in  the  ice  surface,  visible  in  photograph  3.1.  The  remaining  correlations  for  January 
are  more  like  that  of  slab  la.  Figure  3.4  shows  the  computed  autocorrelation  for  the 
slabs  taken  in  March.  This  plot  demonstrates  the  marked  increase  in  the  correlation 
length  for  the  smoother  ice  observed  in  the  March  campaign.  Results  of  the 
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correlation  lengths  and  the  standard  deviation  of  the  surface  height  are  shown  in  table 
3.3.  The  averages  for  both  20  January  and  1 1  March  are  used  in  a  later  section  in 
surface  scattering  models  for  comparison  to  measured  results. 


Table  33. 

Results  of  surface  roughness  calculations 

Date  /  slab  number 

surface  roughness  (mm) 

correlation  length  (cm) 

20  January,  slab  la 

4.068 

1.15 

20  January,  slab  lb 

2.49 

0.29 

20  January,  slab  2a 

1.598 

0.25 

20  January,  slab  2b 

1.427 

0.28 

20  January,  slab  3a 

1.212 

0.37 

20  January,  slab  3b 

2.164 

0.5 

20  January  average 

2.16 

0.47 

1 1  March,  slab  la 

1.209 

1.80 

11  March,  slab  lb 

1.502 

1.72 

1 1  March  average 

1.36 

1.76 
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Photograph  3.1.  Reflected  light  photograph  of  vertical  thick  section  la,  taken  on  20 
Jan.  from  CRREL's  lower  pond  facility.  Grid  scale  measures  1  cm  on  a  side. 


Photograph  3.2.  Reflected  light  photograph  of  verticaJ  thick  section  lb,  taken  on  20 
Jan.  from  CRREL's  lower  pond  facility.  Grid  scale  measures  1  cm  on  a  side. 


Autocorrelation 


Figure  3.1.  Autocorrelation  calculations  from  slab  la  and  lb.  characterizin'' 
the  ice  surface  for  the  January  CRRELEX’90  campaign. 
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Photograph  3.3.  Reflected  light  photograph  of  vertical  thick  section  2a.  taken  on  20 
Jan.  from  CRREL's  lower  pond  facility.  Grid  measures  1  cm  on  a  side. 
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Photograph  3.4.  Reflected  light  photograph  of  vertical  thick  section  2b.  taken  on  20 
Jan.  from  CRREL’s  lower  pond  facility.  Grid  scale  measures  1  cm  on  a  side. 


Autocorrelation 


x'(cm) 

Figure  3.2.  Autocorrelation  calculations  from  slabs  2a  and  2b.  characterizing 
the  ice  surface  for  the  January  portion  of  CRRELEX'90. 
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Photograph  3.5.  Reflected  light  photograph  of  vertical  thick  section  3a.  taken  on  20 
Jan.  from  CRREL's  lower  pond  facility.  Grid  scale  measures  1  cm  on  a  side. 
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Photograph  3.6.  Reflected  light  photograph  of  vertical  thick  section  3b.  taken  on  20 
Jan.  from  CRREL's  lower  pond  facility.  Grid  scale  measures  1  cm  on  a  side. 


Autocorrelation 


Figure  3.3.  Autocorrelation  calculations  from  slabs  3a  and  3b,  characterizing 
the  ice  surface  for  the  January  portion  of  CRRELEX'9(). 
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Photograph  3.7.  Reflected  light  photograph  of  vertical  thick  section  taken  on  1 1 
March  from  CRREL’s  lower  pond  facility.  Note  the  much  smoother  surface  than  the 
January  photographs.  This  is  due  to  melting  and  re-freezing  that  occurred  between 
portions  of  the  CRRELEX'90  campaign. 
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4.  Vector  Calibration  and  Measurements  of  Radar 
Backscatter 


Step-frequency  radars  were  used  to  measure  the  scattering  coefficient  of  saline  ice. 
The  radars  were  calibrated  using  a  complex  vector  calibration  technique  described  in 
this  section.  The  methods  for  obtaining  the  scattering  coefficient  from  a  step- 
frequency  radar  are  described  in  this  section  as  well  as  the  method  for  calculating  the 
backscatter  coefficient  when  applying  the  complex  vector  calibration. 


4.1  Step-frequency  radar  principles 


The  step-frequency  radar  systems  described  in  the  previous  section  were  used  to 
measure  backscatter  from  saline  ice.  The  voltage  scattering  coefficient  from  a  given 
target  can  be  obtained  from  the  step-frequency  radar  measurement  in  the  following 
manner.  The  amplitude  and  phase  of  the  continuous- wave  (CW)  measurement  at 
each  of  the  discrete  frequencies  used  by  the  radar  system  are  stored  on  a  disk.  The 
complex  spectrum  (usually  FFT)  of  these  data  is  related  to  the  voltage  scattering 
coefficient,  sk  of  a  target  by  [lizuka,  et  al.,  1984] 

h  -NS*E0 

hk  z2 


(4.1) 


where  hk  is  the  complex  spectral  component  of  the  measured  data  at  point  k,  N  is  the 
total  number  of  discrete  points  used  over  the  bandwidth,  Eo  is  the  incident  electric 
field  magnitude  and  zfc  is  the  distance  from  the  radar  to  the  scattering  point  in 
question.  The  phase  term  is  simply  an  unwrapping  due  to  the  propagation  of  the 
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incident  and  reflected  wave  and  is  dependent  on  the  specific  frequency  of  the  lc1*1 
measurement  (kfoAz)  as  well  as  the  velocity  (c)  of  the  wave. 

To  obtain  the  power  returned  from  a  target  the  magnitude  of  the  voltage  scattering 
coefficients  are  squared.  This  is  then  used  in  the  narrow-beam  approximation  to  the 
radar  equation  as  the  returned  power. 


4.2  Complex  vector  calibration  technique 

The  radar  return  from  several  calibration  targets  were  used  to  compute  the  vector 
calibration  matrix  as  described  here.  The  targets  included  a  Luneburg  lens,  a  metal 
sphere,  and  a  dihedral  comer  reflector.  The  diplane  was  chosen  so  that  it  could  be 
used  as  both  a  depolarizing  and  a  non-depolarizing  target,  depending  on  its 
orientation.  The  sphere  was  used  for  the  vector  error  calibration  as  a  non¬ 
depolarizing  target.  The  Luneburg  lens  was  used  for  absolute  (scalar)  calibration  of 
the  measurements.  Absolute  calibration  was  also  performed  with  reference  to  the 
sphere  instead  of  the  lens  as  a  verification  of  the  Luneburg  lens  radar  cross  section 
value. 

The  vector  calibration  technique  follows  that  demonstrated  by  [Riegger  and 
Weisbeck,  1989]  for  simple  targets.  Here  the  technique  is  applied  to  the  radar  returns 
from  artificial  sea  ice,  a  complex  distributed  target.  This  vector  calibration  removes 
the  systematic  errors  induced  by  the  radar  system,  including  the  antennas. 
Contributions  to  the  systematic  errors  in  radar  measurements  are  due  to  two  major 
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sources:  frequency  response  of  the  RF  components  of  the  system  and  polarization 
coupling  effects.  The  effects  of  these  error  sources  are  reduced  with  this  technique. 

This  calibration  is  based  on  the  assumption  that  the  measured  voltage  scattering 
coefficient  of  an  object  at  each  polarization  is  the  linear  combination  of  the  true 
voltage  scattering  coefficients  from  the  four  polarizations.  For  example,  the 
measured  voltage  scattering  coefficient  for  VV  polarization,  s™v  is  given  by, 

Sw  =  C10  ■*"  CllSw  +  C12Shh  +  C13Svh  +  C14Shv  (4.2) 

where  s^v,  s^,  Syh,  and  s^v  are  the  true  voltage  scattering  coefficients  for  the  target 
for  each  polarization  combination.  The  nature  of  the  coefficients  Cy  will  be  discussed 
later. 


For  the  set  of  four  linear  polarization  combinations  there  are  three  additional 
equations  similar  to  (4.2).  The  c,q  terms  represent  the  isolation  error,  usually  due  to 
the  anechoic  chamber  in  which  data  are  taken.  If  the  isolation  terms  are  considered  to 
be  negligible  the  following  matrix  equation  results, 


(4.3) 


The  4x4  matrix  is  called  the  calibration  matrix,  C.  The  diagonal  terms  are  due  to 
frequency  response  of  the  components  within  the  system  and  mismatches  within  the 
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system.  These  terms  also  account  for  the  errors  due  to  the  different  physical  paths  for 
the  different  polarizations  in  the  radar  system.  The  off-diagonal  terms  are  due  to 
polarization  coupling  or  polarization  leakage  effects  in  the  transmit  and  receive  paths 
of  the  radar  system.  Leakage  between  polarizations  can  come  from  imperfect 
switches  in  the  RF  section  as  well  as  from  antennas  with  limited  polarization  purity. 
The  amount  of  contribution  to  the  off-diagonal  terms  caused  by  the  switches  is 
determined  by  their  isolation.  The  coefficient  C32,  for  example,  accounts  for  the 
contribution  of  the  HH  component  of  in  the  measurement  of  the  VH  component  of 
(A  This  value  is  roughly  equivalent  to  the  polarization  purity  of  the  antennas  when 
ideal  switches  are  used. 

All  of  these  coefficients  can  be  estimated  by  using  only  three  calibration  targets  if 
several  simplifying  assumptions  are  made.  First,  the  isolation  error  is  ignored  since 
this  is  simply  the  return  with  no  target  In  an  indoor  measurement  facility  the 
chamber  is  primary  the  source  for  these  errors.  In  this  study  the  measurements  were 
made  outdoors  over  a  frozen  saline  pond  and  the  isolation  error  is  negligible.  In 
addition  to  this  assumption,  the  coefficients  C12,  C21,  C34  and  C43  are  assumed  to  be 
negligible.  These  terms,  called  double  leakage  terms,  arise  from  polarization  leakage 
in  both  the  transmit  and  receive  paths  simultaneously.  These  differ  from  the 
remaining  off-diagonal  terms  which  are  due  to  polarization  leakage  either  in  the 
receive  or  the  transmit  path.  The  double  leakage  terms  are  generally  small  in 
comparison  to  the  rest  of  the  matrix  for  antennas  with  modest  polarization  purity.  For 
example,  an  antenna  with  a  polarization  purity  of  20  dB  will  result  in  double  leakage 
terms  that  are  40  dB  down  from  the  dominant  terms,  whereas  the  remaining  off- 
diagonal  terms  are  roughly  20  dB  down. 
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With  these  assumptions  the  remaining  coefficients  can  be  determined  by  performing 
fully  polarimetric  measurements  on  the  three  independent  calibration  targets  .  The 
metal  sphere,  the  diplane  in  a  depolarizing  orientation  (45°  or  135°)  and  the  diplane  in 
a  non-depolarizing  orientation  (0°  or  90°)  are  the  targets  used  to  perform  this 
calibration.  These  targets  can  be  used  to  find  the  remaining  coefficients  and  thus 
quantify  C.  The  sphere  has  a  theoretical  scattering  matrix  given  by,  Ssphere  =  [110 
0  ]T,  where  the  notation  here  for  the  scattering  matrix  is  S  =  [svv,shh,svh,shv]T .  The 

scattering  matrices  for  the  diplane  oriented  vertically  and  horizontally  (0°  and  90°) 
are  given  as,  S  dipianeO  =  [  1  *1  0  0]  ,  and  S  dipiane90  =  [  ~  1  1  0  0].  The  scattering 
matrices  for  the  diplane  in  its  depolarizing  orientations  are,  Sdipiane45  =  [  0  0-1  -1  ]T 
and  Sdipianei35  =  [001  1  ]T,  where  the  angle  is  measured  relative  to  vertical 
direction.  Using  these  known  ("true")  scattering  matrices  and  the  actual  measured 
scattering  matrices,  the  calibration  matrix  C  can  be  obtained.  The  algorithmic  details 
in  computing  C  are  outlined  in  section  4.3. 

This  calibration  matrix  is  then  inverted  so  that  better  estimates  of  the  true  scattering 
coefficients  can  be  obtained  from  the  measured  scattering  coefficients.  The  true 
scattering  coefficients  can  be  obtained  by  a  matrix  multiplication  of  the  inverse  of  the 
scattering  matrix  C  with  the  set  of  measured  scattering  coefficients  from  the  four 
polarization  combinations, 

S°=C"‘  S”.  (4.4) 

Here  Sm  represents  the  measured  scattering  coefficients  and  S°  represents  the  true 
scattering  coefficients  at  the  four  polarization  combinations  as  in  (4.3).  With  this 
vector  calibration  any  fully  polarimetric  measurement  made  by  the  radar  system  can 
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be  used  to  obtain  better  estimates  of  the  true  scattering  coefficients.  The  final  step 
involves  using  thetrue  scattering  coefficients  and  the  radar  equation  for  surface 
targets  to  obtain  a  calibrated  value  for  o°,  the  backscatter  coefficient.  Similarly,  the 
traditional  radar  equation  can  be  used  with  the  corrected  scattering  coefficients  from 
hard  targets  to  obtain  calibrated  radar  cross  sections,  a  for  the  targets.  The  absolute 
calibration  was  made  by  feeding  the  measured  scattering  coefficients  from  either  the 
Luneburg  lens  or  the  sphere  through  the  correction  matrix  procedure  to  obtain  an 
absolute  level.  The  lens  was  used  for  absolute  (scalar)  calibration  because  its 
experimental  radar  cross  section  was  documented  [Bredow,  1989].  Before 
application  to  artificial  sea  ice  returns,  the  vector  calibration  was  applied  to  the  simple 
calibration  targets  in  order  to  validate  the  algorithm  and  reveal  any  limitations. 

4  J  Implementation  of  the  complex  vector  calibration  algorithm 
The  calibration  matrix  C,  which  is  needed  to  characterize  the  radar  system  is 
obtained  by  using  three  calibration  targets.  The  implementation  of  the  algorithm 
begins  with  equation  (4.3).  For  the  sphere  return,  the  components  s°h  ands°v  are 

ideally  zero  so  the  measured  co- polar  terms  can  be  expressed  as 

Sw=cn+Ci2»  (4.5) 

and 

s£,  =c2l+cn.  (4.6) 

When  the  double-leakage  terms  (c  12  and  c2l)  are  neglected,  the  ci  1  and  C22 
coefficients  simply  become  s”  and  sJJ, ,  the  measured  co- polar  scattering  coefficients 

from  the  sphere. 

The  measured  cross-polarized  scattering  coefficients  from  the  sphere  yield  the 
equations, 
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svh  =  C3,  +  C32, 


(4.7) 


and 

Shv  =  C41  +  c42 .  (4.8) 

Likewise,  the  measured  cross-polarized  returns  from  the  diplane  at  0°  orientation 
yield, 

Svh  =  C31  “  C32*  (4.9) 

and 

Sfcv~c4,-c42.  (4.10) 

Equations  (4.7)  and  (4.9)  are  solved  simultaneously  to  obtain  the  coefficients  C31  and 
C32  while  (4.8)  and  (4.10)  are  used  to  obtain  the  coefficients  C41  and  C42.  The 
measurements  from  the  diplane  at  0°  orientation  and  metal  sphere  are  used  to 
calculate  the  left  half  of  the  calibration  matrix.  If  the  double-leakage  terms  cannot  be 
neglected  they  may  be  obtained  by  solving  the  set  of  simultaneous  equations  obtained 
from  the  co- polar  measurements  from  both  the  sphere  and  diplane  instead  of  just  the 
sphere  return  as  outlined  above. 

The  diplane  at  45°  orientation  is  used  in  the  calculation  of  the  C33  and  C44  returns. 

Since  the  true  co-polar  scattering  coefficients  for  the  diplane  in  this  orientation  are 
zero  the  following  equations  relate  the  measured  returns  to  the  calibration 
coefficients. 


;vh  -  (C33  +C3<)» 

(4.11) 

'tav  =  ~(C43  +C44)' 

(4.12) 

41 


If  the  double-leakage  terms  (C34,  C43)  are  once-again  neglected  the  coefficients  C33 
and  C44  are  simply  the  measured  cross-polarized  returns  from  the  diplane  at  45° 
orientation. 

The  remaining  terms  to  be  computed  are  the  coefficients  ci3,  ci4,  C23,  and  C24. 

These  are  calculated  with  the  measured  co- polarized  return  when  the  diplane  is  at  45°. 
The  assumption  required  for  this  is  that  the  coefficients  c  13  and  c  14  are  equivalent,  as 
are  the  coefficients  C23  and  C24-  This  yields 

Cn=c14=*s“,  (4.13) 

and 

c23=cM=^s”  (4.14) 

as  the  estimates  for  these  coefficients.  While  the  assumption  that  these  pairs  of 
coefficients  are  equivalent  seems  rather  crude  on  the  surface,  it  holds  if  the  cross- 
polarized  scattering  coefficients  to  be  measured  are  the  same. 

Figure  4.1  shows  the  steps  involved  in  applying  this  calibration  technique  to 
determine  the  scattering  coefficients  from  a  distributed  targets  such  as  a  saline  ice 
sheet.  A  parallel  procedure  can  be  used  to  "re-calibrate"  the  radar  cross  sections  of 
the  calibrations  themselves  in  order  to  verify  the  effectiveness  of  the  technique  as 
well  as  the  absolute  level  of  the  calibration.  The  schematic  of  this  "re-calibration" 
procedure  is  shown  in  Figure  4.2.  This  was  used  to  cross-calibrate  the  radar  with  the 
sphere  and  lens  measurements  to  obtain  the  most  accurate  absolute  calibration 
possible. 
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Figure  4.1.  Details  of  the  calibration  procedure. 


Correct  Absolute 
a  Reference 


Figure  4.2.  Measure  of  radar  cross  sections  of  hard  targets 
using  vector  calibration  procedure 
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4.4  Narrow-beam  approximation  to  o° 


For  the  usual  scalar  calibration,  the  radar  equation  is  modified  in  order  to  account  for 
a  distributed  target.  When  the  target  in  question  is  a  distributed  one  the  radar 
equation  can  be  written  as 


P 


pq 


(4.15) 


If  the  antenna  gain  pattern  is  assumed  to  be  a  narrow  beam  then  the  gain  is 

approximated  as  constant  over  the  area  illuminated.  The  range  as  well  as  the 

variation  also  are  assumed  to  vary  little  across  the  beamwidth  yielding, 

q°  _  Prge»l^  P<flc«l 

Pc^AP* 


(4.16) 


Here  Pr  is  the  power  returned  from  the  distributed  target  and  PCal  is  the  power 
returned  from  a  calibration  target  of  known  radar  cross  section,  aCal-  The  parameters 
Pdlcal  and  Pdl  are  the  delay  line  readings  taken  at  the  time  of  calibration  and  the  time 
of  the  field  measurement,  respectively.  The  range  to  the  distributed  target  is  R,  A  is 
the  area  illuminated  by  the  antenna,  and  Real  is  the  range  to  the  calibration  target. 


This  equation  for  the  calculation  for  the  scattering  coefficient  is  often  called  the 
narrow-beam  approximation.  The  effect  of  the  antenna  pattern  and  the  variation  of  o° 
across  the  illuminated  area  is  taken  into  account  in  a  post-processing  step  using  the 
technique  described  by  [Wang  and  Gogineni,  1991]. 

The  narrow-beam  approximation  also  assumes  that  the  polarization  purity  of  the 
antennas  is  infinite,  that  is  there  is  no  horizontal  component  of  the  electric  field 
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transmitted  (or  received)  when  the  system  is  switched  to  vertical  polarizations.  If  the 
polarization  purity  is  finite,  as  in  any  real  antenna  system,  the  integral  equation  for  the 
scattering  coefficient  becomes  more  complicated.  The  effect  of  the  antenna 
limitations  is  demonstrated  by  the  following.  Let  Gvv  be  the  gain  of  the  vertical 
component  of  the  antenna  system  in  vertical  mode,  Gvh  be  the  gain  of  the  vertical 
component  of  the  antenna  system  in  horizontal  mode,  Ghv  be  the  gain  of  the 
horizontal  component  of  the  antenna  system  in  vertical  mode,  and  Ghh  be  the  gain  of 
the  horizontal  component  of  the  antenna  system  in  horizontal  mode. 

Then  the  measured  scattering  coefficient  in  terms  of  the  true  scattering  coefficients 
can  be  expressed  as. 


p  12a 

P  -  rfK  2  0  q2  0  +G2  0  +(J2  0  j 

PP  (4jj)3R4  1  »  PP  «  «  M  M  ^qp^qpj- 


(4.17) 


For  example  the  returned  power  for  the  W  component  is  given  as 


.1 2*hi 


p  -  rr.2 


(4ti)  R 


L[G^<C  +GS.OS,  +G>1  +G*,cCl. 


(4.18) 


This  is  the  scalar  approximation  to  the  vector  model  assumed  in  the  complex  vector 
calibration  approach,  where  the  measured  <A  is  expressed  as  the  linear  combination 
of  the  true  scattering  coefficients.  If  we  used  the  voltage  scattering  coefficients, 
which  are  related  to  O0  by  s  <*  ,  then  the  matrix  formulation  of  the  calibration 

matrix  results.  The  coefficients  of  the  matrix  are  complex  in  order  to  account  for  both 
the  magnitude  and  phase  of  the  field  quantities.  As  mentioned  before,  the  complex 
voltage  scattering  coefficient  must  be  squared  as  a  measure  of  the  power  returned 
from  the  target  to  be  used  in  equation  4.16  to  obtain  a  measure  of  <A 


4.5  Beamwidth  correction 

The  use  of  wide-beam  antennas  results  in  a  second  source  of  systematic  error  if  the 
narrow-beam  approximation  (equation  4.16)  is  used  in  calculating  o°.  Since  the 
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measurements  were  taken  with  relatively  wide-beam  antennas  and  the  calculations 
used  for  o°  relied  on  the  narrow-beam  approximation,  a  wide-beam  correction  had  to 
be  performed  on  the  data.  In  [Wang  and  Gogineni,  1991]  a  numerical  procedure  for 
recovering  the  backscatter  coefficient  from  wide-beam  measurements  was  presented. 
This  technique  utilizes  a  two-part  linear  model  fit  to  the  experimental  data  as  well  as 
the  assumption  of  a  Gaussian  antenna  beam.  These  assumptions  are  used  to  solve  the 
integral  form  of  the  radar  equation  for  surface  scattering.  The  difference  between  the 
assumed  model  and  the  results  of  integration  are  used  to  correct  the  actual  measured 
data.  Because  of  this,  the  technique  is  relatively  insensitive  to  the  model  chosen. 

This  procedure  was  applied  to  the  CRRELEXIX)  lower  pond  data  set  and  the  results 
are  shown  in  the  following  section. 
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5.  CRRELEX'90  Results 


The  calibration  was  applied  to  several  targets  to  verify  its  effectiveness.  In  addition 
to  the  targets  used  to  generate  the  vector  error  calibration  matrix,  a  Luneburg  lens  and 
the  diplane  in  a  second  orientation  were  used.  The  vector  calibration  was  then 
applied  to  the  artificial  sea  ice  returns  to  obtain  high  quality  radar  backscatter 
measurements.  The  radar  cross  sections  of  the  hard  targets  as  well  as  the  o°  of  the  ice 
were  measured  with  the  normal  scalar  calibration  as  well  for  comparison. 

5.1  Calibration  target  results 

The  vector  calibration  algorithm  was  applied  to  the  calibration  targets  in  order  to 
validate  the  technique.  The  theoretical  radar  cross  sections  for  the  calibration  targets 
are  shown  in  table  5.1  for  both  the  targets  at  both  Ku  band  and  C  band. 


Table  5.1.  Theoretical  Radar  Cross  Sections  for  Calibration  Targets 

Ku-band  a  (dBm^) 

C-band  o  (dBm^) 

Calibration  Target 

1 

Luneburg  Lens 

14.7 

-oo 

6.7 

-oo 

Metal  Sphere 

-13.9 

-oo 

-13.9 

-oo 

Diplane  (0°  or  90°) 

23.3 

-oo 

16.2 

-oo 

Diplane  (45°  or  135°) 

-oo 

23.3 

-OO 

16.2 

The  results  for  the  measurements  of  the  calibration  targets  at  Ku  band  are  shown  in 
table  5.2  and  the  results  from  the  C-band  measurements  are  in  table  5.3.  Here,  the 
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calculation  of  o  is  shown  for  the  case  when  only  simple  scalar  calibration  is  used  and 
when  the  complex  vector  calibration  is  applied. 


Table  5.2.  Ku-band  Radar  Cross  sections  in  dBm2 


Calibration 

Scalar  Calibration 

Vector  Calibration 

Target 

°w 

°w 

°vh 

°hv 

sphere 

-12.97 

-10.9 

-25.3 

-22.8 

-13.9 

-13.9 

-51.9 

-52.1 

Diplane  0° 

17.5 

18.6 

-2.82 

-1.82 

16.6 

15.6 

-22.0 

-22.2 

Diplane  45° 

4.58 

2.52 

18.7 

20.4 

-322 

-312 

17.9 

17.9 

Diplane  90° 

23.1 

24.5 

-0.68 

-0.51 

22.1 

21.5 

-14.4 

-18.1 

Lens 

14.9 

12.8 

-5.77 

-9.92 

13.9 

9.47 

-15.6 

-14.5 

Table  5 .3.  C-band  Radar  Cross  sections  in  dBm2 


Calibration 

Scalar  Calibration 

Vector  Calibration 

Target 

®w 

°hh 

°vh 

^hv 

Oyv 

°vh 

°hv 

sphere 

-20.4 

-14.5 

-39.2 

-37.3 

-20.4 

-15.9 

-63.7 

-68.7 

Diplane  0° 

6.90 

5.81 

-9.76 

-14.0 

6.93 

4.36 

-36.9 

-41.2 

Diplane  45° 

-15.0 

-8.07 

7.76 

8.09 

-341 

-338 

7.77 

7.34 

Diplane  90° 

9.15 

10.1 

-8.15 

-8.57 

9.10 

8.87 

-24.0 

-39.3 

Lens 

6.70 

8.14 

-11.9 

-24.6 

6.70 

6.70 

-39.3 

-35.8 

The  ratio  of  co-polarized  o  to  cross-polarized  o  is  increased  drastically  for  the  sphere, 
lens  and  diplane  at  0°  and  diplane  at  90°.  These  are  ideally  non-depolarizing  targets. 
These  results  demonstrate  that  the  vector  error  correction  technique  is  able  to  remove 


systematic  errors  and  obtain  a  better  estimate  of  the  true  res  of  these  targets.  For  the 
diplane  at  45°,  an  ideally  completely  depolarizing  target,  the  cross-polarized  to  co¬ 


polarized  res  ratio  is  also  improved  dramatically.  This  is  illustrated  by  figures  5.1  and 
5.2,  which  show  the  polarization  ratio  for  the  targets  measured.  The  polarization  ratio 


for  ideally  non-depolarizing  targets  is  defined  as, 

p  =  °w  +  ghh 


^vh  °hv 


(5.1) 


and  for  ideally  depolarizing  targets  as, 


p  —  ^vh  ^hv 

<>w +<*,*' 


(5.2) 


Here,  Ojj  is  the  radar  cross  section  for  ii  polarization  combination  and  i  or  j  denotes  v 
or  h,  representing  vertical  or  horizontal  polarization.  The  polarization  ratio  for  each 
of  the  hard  targets  when  the  vector  calibration  is  applied,  compared  to  the  scalar 
calibration  is  shown  in  figure  5.1  for  the  Ku-band  case.  The  C-band  results  are  shown 
in  figure  5.2.  These  results  demonstrate  the  ability  of  the  vector  calibration  to 
effectively  increase  the  dynamic  range  of  the  measurements  performed  by  the  radar 
systems. 
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Ku-band  Calibration 
Polarization  Ratios 


Diplans  (angle  of  orientation) 


Figure  5.1.  Polarization  ratio  comparison  for  Ku-band  targets  between  scalar 
and  vector  calibration.  The  polarization  ratio  for  the  diplane  at  45°  is  limited 
by  the  dynamic  range  used  in  the  computation  since  the  calculation  depends 
on  this  return. 


C-band  Calibration 
Polarization  Ratios 


Figure  5.2.  Polarization  ratio  comparison  for  C-band  targets  between  scalar 
and  vector  calibration. 
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5.2  Polarimetric  responses 

Polarimetric  response  of  the  calibration  targets  were  derived  from  the  fully 
polarimetric  data  before  and  after  applying  the  vector  calibration  to  the  data.  The 
polarimetric  response  (sometimes  called  polarization  signature)  gives  the  scattering 
coefficient  for  an  arbitrary  incident  polarization.  The  co-polarized  response  gives  the 
scattering  coefficient  for  the  same  polarization  as  would  be  transmitted  as  a  function 
of  the  transmitted  polarization.  The  cross-polarized  response  gives  the  scattering 
coefficient  of  the  orthogonal  polarization  as  a  function  of  the  incident  polarization. 
The  polarization  is  defined  by  the  ellipticity  and  the  orientation  of  the  major  axis  with 
respect  to  vertical.  The  polarization  response  can  be  obtained  from  the  measurement 
of  the  (voltage)  scattering  coefficient  at  the  four  linear  polarizations:  VV,  HH,  VH, 
and  HV  [Riegger,  1989]. 


The  scattering  matrix  for  any  pair  of  similarly  related  orthogonal  polarizations  can  be 


obtained  from  the  linear  ones  by  the  relation, 

[s-]  =  [Q(a2,'F2)]-,[SKQ(a1,'F,)], 

where. 


and  SA  is  a  new  coordinate  basis, 


Q(ai,\|/i)  is  an  operator  associated  with  the  incident  field  and  Q(oc2,\|/2)  is  the 


(5.3) 


(5.4) 


(5.5) 


operator  associated  with  the  scattered  field.  These  are  given  in  terms  of  the  angles  a 


and  \j f  as, 
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[" cos  y  cos  a  -  j  sin  y  sin  a  -  sin  y  cos  a  +  j  cos  y  sin  a' 
^  sinycosa  + jcosysina  cosycosa  + jsinysina 


(5.6) 


The  angle  oc  corresponds  to  the  axial  ratio  of  the  polarization  ellipse  describing  either 
the  incident  (with  subscript  1)  or  scattered  (subscript  2)  polarization.  For  linear 
polarizations  a=0°,  and  for  circular  polarizations  a=±45°.  The  angle  y  is  the  amount 
of  rotation  in  the  plane  perpendicular  to  the  direction  of  propagation  that  the  wave 
polarization  vector  must  move  to  align  with  the  vertical  direction  (y=0°).  The  new 
scattering  coefficient  s^  plotted  versus  a  and  y  define  what  is  commonly  referred  to 
as  the  co-polarized  polarization  response.  This  is  a  plot  of  the  co-polarized  scattering 
coefficient  for  an  arbitrary  incident  polarization  vector  and  corresponding  scattering 
polarization.  Similarly,  the  plot  of  Si2  versus  these  angles  is  the  cross-polarized 
response.  This  is  a  scattering  coefficient  plot  versus  an  arbitrary  incident  polarization 
with  the  corresponding  orthogonally-polarized  scattered  field. 


5.3  Polarimetric  responses  for  the  calibration  targets 

The  polarization  responses  for  the  calibration  targets  measured  during  CRRELEX'90 
were  derived  from  the  linear  polarization  measurements.  The  ideal  polarization 
responses  for  the  calibration  targets  used  are  shown  in  figures  5.3,  5.4  and  5.5  [Ulaby, 
1990] .  These  are  used  to  judge  the  quality  of  the  vector  calibration  as  it  applies  to 
the  calibration  targets.  The  diplane  polarization  response  demonstrates  the  value  of 
the  orientation  variable  in  the  signature.  When  the  diplane  is  rotated  by  45°  the 
polarization  response  shows  a  shift  of  45°  along  the  orientation  axis  as  shown  in 
figures  5.4  and  5.5.  This  is  useful  in  radar  classification  of  targets  when  the 
orientation  is  unknown. 
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The  polarization  responses  for  the  calibration  targets  were  computed  before  and  after 
applying  the  vector  calibration.  These  are  shown  in  the  following  series  of  figures  for 
C  band,  followed  by  the  same  for  Ku  band.  For  all  figures,  the  co-polarized 
polarization  response  is  on  the  left  and  the  cross-polarized  response  on  the  right. 
These  demonstrate  the  improved  responses  for  the  calibration  targets  when  the  vector 
calibration  is  used. 


Figure  5.3.  Ideal  polarization  response  for  a  perfectly-conducting  metal  sphere.  The 
co-polarized  polarization  response  is  on  the  left  and  the  cross-polarized  response  is  on 

the  right. 


CO 


CR 


Figure  5.4.  Ideal  polarization  response  for  a  dihedral  comer  reflector  oriented  either 
vertically  or  horizontally.  The  co- polarized  polarization  response  is  on  the  left  and 
the  cross-polarized  response  is  on  the  right. 


Figure  5.5.  Ideal  polarization  response  for  a  dihedral  comer  reflector  oriented  45° 
from  vertical.  The  co- polarized  polarization  response  is  on  the  left  and  the  cross- 
polarized  response  is  on  the  right.  Note  the  shift  of  45°  along  the  orientation  axis 

from  figure  5.4. 
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Figure  5.6  C-band  polarimetric  response  for  the  metal  sphere  before  the  complex 
vector  calibration  was  applied.  The  co-polarized  response  is  on  the  left  and  the  cross- 

polarized  response  is  on  the  right. 


Figure  5.7  C-band  polarimetric  response  for  the  metal  sphere  after  the  complex 
vector  calibration  was  applied.  The  co-polarized  response  is  on  the  left  and  the  cross- 

polarized  response  is  on  the  right. 
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Figure  5.8  C-band  polarimetric  response  for  dihedral  comer  reflector  at  0° 
orientation  before  the  complex  vector  calibration  was  applied. 


Figure  5.9  C-band  polarimetric  response  for  dihedral  comer  reflector  at  0° 
orientation  after  the  complex  vector  calibration  was  applied. 
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Figure  5.10  C-band  polarimetric  response  for  dihedral  comer  reflector  at  45° 
orientation  before  the  complex  vector  calibration  was  applied.  The  co- polarized 
response  is  on  the  left  and  the  cross-polarized  response  is  on  the  right. 


Figure  5.11  C  -band  polarimetric  response  for  dihedral  comer  reflector  at  45° 
orientation  after  the  complex  vector  calibration  was  applied.  The  co-polarized 
response  is  on  the  left  and  the  cross-polarized  response  is  on  the  right. 
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Figure  5.12  C-band  polarimetric  response  for  dihedral  comer  reflector  at  90° 
orientation  before  the  complex  vector  calibration  was  applied.  The  co-polarized 
response  is  on  the  left  and  the  cross-polarized  response  is  on  the  right. 
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Figure  5.13  C-band  polarimetric  response  for  dihedral  comer  reflector  at  90° 
orientation  after  the  complex  vector  calibration  was  applied.  The  co-polarized 
response  is  on  the  left  and  the  cross-polarized  response  is  on  the  right. 
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Figure  5.14  Ku-band  polarimetric  response  for  the  metal  sphere  before  the  complex 

vector  calibration  was  applied. 
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Figure  5.15  Ku-band  polarimetric  response  for  the  metal  sphere  after  the  complex 

vector  calibration  was  applied. 
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Figure  5.16  Ku-band  polarimetric  response  for  the  diplane  at  0°  orientation  before 
the  complex  vector  calibration  was  applied. 


Figure  5.17  Ku-band  polarimetric  response  for  the  diplane  at  0°  orientation  after 
applying  the  complex  vector  calibration. 
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Figure  5.18  Ku-band  polarimetric  response  for  the  diplane  at  45°  orientation  before 
applying  the  complex  vector  calibration. 


Figure  5.19  Ku-band  polarimetric  response  for  the  diplane  at  45°  orientation  after 
applying  the  complex  vector  calibration. 
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Figure  5.20  Ku-band  polarimetric  response  for  the  diplane  at  90°  orientation  before 
applying  the  complex  vector  calibration. 
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Figure  5.21  Ku-band  polarimetric  response  for  the  diplane  at  90°  orientation  after 
applying  the  complex  vector  calibration. 
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Note  the  dramatic  improvement  in  the  polarization  responses  for  the  diplane  at  90° 
orientation  and  the  sphere  at  C  band.  The  other  polarization  responses  demonstrate 
some  improvement,  which  is  evident  in  the  raising  and  sharpening  of  the  peaks  and 
valleys  in  the  responses. 


5.4  Simulated  sea  ice  results 

Once  the  effectiveness  of  the  method  was  verified  the  calibration  was  applied  to 
complex  returns  from  the  artificial  sea  ice  measurements.  Two  resulting  backscatter 
coefficient  (o°)  measurements  for  data  taken  during  the  March  experiment  are  shown 
in  figures  5.22  and  5.23  for  the  C-band  and  Ku-band  systems,  respectively.  For 
comparison,  o°  measurements  using  the  scalar  calibration  are  also  shown  in  figures 
5.24  and  5.25.  These  demonstrate  the  improvement  in  dynamic  range  of  the 
measurements  when  the  vector  calibration  is  applied.  The  ratio  of  co- polar  to  cross- 
polar  a0  is  much  higher  for  the  vector-calibrated  data  than  for  the  scalar-calibrated 
data.  These  results  are  corrected  for  errors  introduced  by  wide-beam  antennas  as 
described  in  section  4.5. 

These  results  are  from  the  early  March  portion  of  the  CRRELEX'90  campaign.  These 
curves  are  characteristic  of  those  from  slightly  rough  surfaces.  Both  the  vector 
calibration  and  the  scalar  calibration  were  applied  to  all  measurements  made  during 
the  experiment  at  CRREL’s  lower  pond  facility.  The  resulting  o°  vs.  angle  of 
incidence  for  both  the  scalar  and  vector  calibration  schemes  are  shown  in  appendices 
A  and  B  before  applying  the  wide-beam  antenna  correction.  Appendix  A  shows  the 
scalar  calibrated  results  for  Ku  band  and  C  band.  Appendix  B  shows  the  Ku-band 
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and  C-band  results  after  applying  the  vector  calibration.  The  resulting  o°  vs.  angle  of 
incidence  plots  after  the  application  of  the  wide-beam  antenna  correction  are  shown  in 
Appendix  C  for  scalar  calibration  and  in  Appendix  D  for  the  vector  calibration. 

These  backscatter  measurements  show  a  strong  angular  trend,  characteristic  of 
scattering  from  smootn  surfaces.  The  results,  after  careful  calibration,  are  compared 
with  surface  scattering  theory  in  the  next  section.  This  demonstrates  that  the  primary 
backscatter  mechanism  from  the  saline  ice  during  CRRELEX'90  was  surface  scatter. 
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CRRELEX 
C  Band  Mar.  8,  '90 
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Figure  5.22  Backscatter  from  saline  ice  at  C  band  using  the  complex  vector 
calibration.  These  measurements  show  the  improvement  in  the  dynamic  range  of  the 

measurement 
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Ku  Band  Mar.  5,  '90 
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Figure  5.23  Backscatter  from  saline  ice  at  Ku  band  using  the  complex  vector 
calibration.  These  measurements  show  the  improvement  in  the  dynamic  range  of  the 

measurement 
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CRRELEX 
C  Band  Mar.  8,  '90 


Incidence  Angle 

Figure  5.24  Backscatter  from  saline  ice  at  C  band  using  scalar  calibration.  These 
measurements  show  the  effect  of  polarization  leakage  in  the  cross-polarized  return. 

CRRELEX 
Ku  Band  Mar.  5, '90 
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Incidence  Angle 


Figure  5.25.  Backscatter  from  saline  ice  at  Ku-band  using  scalar  calibration.  These 
measurements  show  the  effect  of  polarization  leakage  in  the  cross-polarized  return. 
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6.  Comparison  of  Backscatter  Measurements  From 

Simulated  Sea  Ice  to  Rough  Surface  Scattering  Theory 

The  microwave  backscatter  measurements  performed  at  CRREL  were  fit  to  simple 
scattering  models  to  verify  the  quality  of  the  data  set  and  to  demonstrate  that  the 
primary  backscatter  mechanism  is  surface  scattering.  The  first  of  these  models  is  the 
small  perturbation  model  for  backscatter  from  a  non-periodic  random  rough  surface. 

A  random  rough  surface  is  characterized  by  the  dielectric  constant  of  the  medium,  the 
surface  correlation  function,  and  the  standard  deviation  of  the  surface  height  Oh,  also 
called  the  rms  surface  height.  The  correlation  length,  Lc  is  a  commonly  used 
parameter  extracted  from  the  correlation  function.  The  correlation  lengths  and 
surface  roughness  were  computed  for  the  CRREL  experiment  as  shown  in  section  4. 
A  short  table  of  the  average  surface  parameters  is  shown  in  table  6.1 


Table  6.1  CRRELEX'90  sea  ice  surface  parameters 

Date 

Oh  (mm) 

Lcorr  (mm) 

20  Jan 

2.16 

4.7 

11  Mar 

1.36 

17.6 

6.1  Random  rough  surface  scattering 

The  scattering  from  any  surface  consists  of  two  components:  a  coherent  and  a  non¬ 
coherent  component.  The  total  backscatter  is  the  sum  of  the  incoherent  and  coherent 
components. 
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Several  theories  to  solve  the  incoherent  scattering  from  rough  surfaces  exist  in  the 
literature,  including  integral  equation  methods  [Fung  and  Pan,  1987],  the  Kirchoff 
model  [Rice,  1951]  and  [Wu  and  Fung,  1972],  and  the  small  perturbation  method 
[Ulaby,  et  al.,  1982].  The  most  commonly  used  are  the  Kirchoff  model  and  the  small 
perturbation  method.  The  Kirchoff  model  uses  the  scalar  approximation  to  the  wave 
equation  in  order  to  find  the  scattered  field  intensities.  This  is  considered  a  physical 
optics  (PO)  approximation  and  applies  to  high  frequency  cases,  when  then  surface 
parameters  are  large  with  respect  to  wavelength.  The  small  perturbation  method 
assumes  that  the  surface  height  variations  are  small  with  respect  to  wavelength.  The 
boundaries  of  limitation  of  the  small  perturbation  method  are  shown  in  table  6.2  and 
the  limits  for  the  Kirchoff  approach  are  shown  in  table  6.3  [Ulaby,  et  al.,  1986].  Also 
in  these  tables  are  the  calculated  parameters  using  the  CRRELEX'90  measurements 
for  both  frequencies  and  both  time  periods  as  an  aid  in  determining  which  scattering 
model  may  be  appropriate. 


Table  6.2  Limitations  of  small  perturbation  metl 

hod. 

Parameter 

Lc 

k*Lc 

k*Oh 

Limits 

<0.212 

<3 

<0.3 

n 

0.46 

0.52 

0.24 

C-band  March 

0.08 

1.95 

0.15 

1— 

0.46 

1.32 

0.61 

Ku-band  March 

0.08 

4.94 

0.38 
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Table  6.3  Limitations  of  Kirchoff  model 

Parameter 

k*Lc 

JL 

ohX 

Limits 

>6 

>2.7 

C-band  January  values 

0.52 

0.21 

C-band  March  values 

1.95 

0.52 

1.32 

0.5 

Ku-band  March  values 

4.94 

10.2 

The  coherent  component  of  the  backscatter  is  due  to  the  direct  reflection  at  the  surface 
and  dominates  near  nadir.  Also  called  the  specular  component,  it  can  be  described  as 
[Ulaby,  et  al.,  1982] 

oc  =  (aof2  5(q  x  )5(qy ) ,  (6.1) 

where  kj  is  the  wave  number  in  the  first  medium,  usually  free  space,  and  ao  is 

ac  =  -2Rh0(cos8)  (6.2) 


for  HH  polarization  and 

a  „  =  -2RvO(cos0)  (6.3) 

for  VV  polarization.  Here  RhO  is  the  perpendicular  (horizontal)  polarization 
reflection  coefficient  and  RvO  is  the  parallel  (vertical)  polarization  reflection 
coefficient.  Due  to  the  dirac  delta  functions,  8(qx)  and  8(qy),  the  coherent  component 
exists  only  for  0°  incidence  angle  in  the  ideal  case.  When  using  finite-beam  antennas 
a  first-order  solution  to  the  contribution  due  to  the  coherent  component  of  the 
scattering  is,  [Fung  and  Eom,  1983] 

^s^^exp^kXJexK-®2/82).  (6.4) 
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where  p  is  the  polarization,  0  is  the  incidence  angle,  Tp  is  the  reflection  coefficient  for 
polarization  p  (either  H  or  V),  and  B  is  given  by 

B!=^+(P/2)I  <( 

The  one-sided  beamwidth  of  the  antennas  is  P  and  the  range  to  the  scatterer  is  Ro* 

The  coherent  component  of  the  scattering  must  be  considered  when  the  rms  surface 
height  is  smaller  than  A./4  [Ulaby,  et  al.,  1986].  The  rms  surface  roughness  is  smaller 
than  A/10  for  Ku  band  and  A/25  for  C  band  in  the  highest  rms  roughness  case  during 
CRRELEX'90.  Because  of  this,  the  coherent  component  of  the  scattering  must  be 
considered. 

6.2  Small  perturbation  method 

The  coherent  component  must  be  combined  with  the  non-coherent  scattering 
component  to  fully  describe  the  backscatter  from  a  slightly  rough  surface.  The 
models  considered  in  this  study  for  describing  the  incoherent  scattering  component 
are  reviewed  in  this  section. 

The  small  perturbation  method  for  determining  the  backscatter  from  a  slightly  rough 
random  surface  is  fully  developed  in  [Ulaby,  1982]  in  general.  The  primary 
assumptions  in  the  small  perturbation  method  are  that  the  surface  standard  deviation, 
oh  is  small  with  respect  to  wavelength  and  the  surface  slope  is  also  small.  These 
conditions  are  shown  in  table  6.2  with  approximate  limits  as  imposed  through 
experimentation  [Ulaby,  1982]  and  [Chen,  1988].  These  assumptions  allow  the  fields 
to  be  approximated  by  the  first  few  orders  of  a  Taylor's  series  expansion.  The  first- 
order  solution  results  are  described  below. 
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The  result  for  backscatter  and  non-magnetic  surfaces  (|ir=l)  reduces  to 

a'n  =  8k4a2cos40jaM|2W(2ksin0,O), 

where  pq  denotes  the  polarization  combination  (VV,  HH,  VH,  or  HV),  k  is  the  wave 
number  in  medium  1,  and  W  is  the  Fourier  Transform  of  the  surface  correlation 
function.  The  parameters  otpq  are  given  by, 

£,"1 

ahh  -  T  / 

cos0  +  ye,  -  sin  0 1 

.  sin20-e.(l  +  sin20) 


rcos  0  +  -Je, 


.  ,  -|2  * 

r  -sin20J 


ahv  =avh  =0. 


This  is  the  first-order  approximation  for  the  cross-polarized  return.  No  cross- 
polarized  component  due  to  surface  scatter  alone  arises  in  small  perturbation  theory 
until  higher  orders  are  considered  [Valenzuela,  1967].  The  technique  for  calculating 
the  theoretical  cross-polarized  return  is  shown  in  section  6.3. 


Two  common  assumptions  for  the  correlation  function  of  random  surfaces  have  been 

used  to  compare  measurements  with  theory.  These  are  the  Gaussian  correlation 

function  and  the  exponential  correlation  function.  It  has  been  shown  through 

experiment  [Bredow,  1989]  that  the  exponential  correlation  function  assumption  fits 

well  to  measured  data  from  saline  ice  grown  at  CRREL.  Using  the  exponential 

correlation  assumption,  equation  6.6  simplifies  to, 

8(kah  )2  (kLc  )2  cos4  0|a 


[l  +  (2kLc  sin©)2] 


This  can  be  easily  computed  for  the  set  of  surface  parameters  as  a  function  of 
incidence  angle  for  VV  and  HH  polarization.  This  is  performed  in  section  6.5  for 
comparison  with  the  CRRELEX'90  measurements.  The  cross-polarized  component  is 
zero  under  the  first-order  approximation.  The  following  section  describes  the  method 
by  which  the  second  order  approximation  to  the  cross-polarized  return  is  obtained. 


63  Cross-polarization  calculations  with  second-order  small  perturbation  method 
The  cross-polarized  backscatter  from  random  rough  surface  comes  from  using  a 
second  order  approximation  of  the  perturbation  series  expansion  [Ulaby,  et  al.,  1986, 
pp.  952].  When  the  second  order  is  taken  into  account,  the  cross-polarized 
backscatter  coefficient  is. 


_o  _  _o  _ 

CTvh  ~  ®hv  — 


-  •*  2  2 
“  ‘  uv 


7tk4Oh  COS2  0 


(e.-lXR.-Rt)f 

2 


(6.8) 


•  f  f  ^-^rW(u-ksin0,v)W(u  +  ksin8,v)dudv 

iiW 


where 

D0  =  er^ko  -  (u2  +  v2)  +  <^erko  -  (u2  +  v2) , 
and  0  is  the  angle  of  incidence.  W(u,v)  is  the  roughness  spectrum  and  is  related  to  the 
surface  correlation  p(a,|3)  by, 

W(u,v)  =  ^  J  J p(a,(3)e(_jau_jPv)dcxdp.  (6.9) 


If  the  integral  is  transformed  into  cylindrical  coordinates  the  final  form  of  the 

solution,  suitable  for  calculation  on  a  computer,  is  given  as: 

R5cos2<)>sin2<t)W(R+)W(R~) 


0 

The  coefficient  term  Tvh  can  be  expressed  as 


<=<v=2Tck4o4cos40T¥h. 


“2* 

If 


(6.10) 
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(6.11) 


and 

R*  =  -^R2  +  koSin20±2koRsin0cos<|>.  (6.12) 

If  the  correlation  function  is  assumed  to  be  exponential,  the  isotropic  roughness 
spectrum  can  be  expressed  as, 

L2 

W(R)  = - =* - 5- .  (6. 1 3) 

[l  +  (LcR)2p 

The  integration  in  6.10  must  generally  be  carried  out  numerically.  The  integrand 
tends  to  be  well-behaved  for  the  range  of  allowable  values  for  oh  an  Lc  allowed  in  the 
small  perturbation  method. 


6.4  Comparison  of  measurements  with  theory 

Comparison  of  the  measured  cross-polarized  response  with  the  theoretical  value 
computed  from  the  theory  presented  here  is  made  in  this  section.  Using  the  small 
perturbation  method  and  the  surface  parameters  computed  for  CRRELEX'90,  the 
theoretical  backscatter  coefficient  were  calculated.  These  computations  were 
compared  with  the  measurements.  Figure  6.1  shows  the  result  of  the  small 
perturbation  method  plotted  along  with  the  measured  Ku-band  co-polarized  response 
and  the  result  of  the  method  of  6.3  with  the  cross- polarized  response.  The  backscatter 
measurements  were  taken  on  the  same  day  as  the  sample  from  which  the  surface 
roughness  parameters  were  obtained.  Since  the  cross-polarized  response  was 
observed  to  be  overestimated  in  this  case,  the  correlation  length  was  adjusted  to 
achieve  a  better  fit  to  the  measure  data.  This  parametric  fit  is  shown  in  figure  6.2  for 
a  correlation  length  of  3.48  cm.  The  correlation  length  computation  from  the  thick 
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section  photographs  may  be  a  poor  estimate  of  the  actual  correlation  length  due  to  the 
small  sample  length  in  the  photographs.  In  slab  la,  for  example,  (section  4,  figure 
4.1)  the  correlation  length  calculation  is  twice  that  of  slab  lb  although  both  were 
taken  at  the  same  time. 


Incidence  Angle 

Figure  6.1  Comparison  of  small  perturbation  method  with  measured 
backscatter  at  Ku  band  from  CRRELEX'90.  Theory  for  the  cross-polarized 
response  is  fie  integral  equation  method.  All  theoretical  curves  are  drawn  as 
lines,  with  the  measurements  as  various  symbols.  This  demonstrates  that  the 
primary  backscatter  mechanism  is  surface  scattering.  Values  used  in  model 
are  Oh =1.36  mm,  Lc=17.6  mm  and  £f=3.75-j0.2. 
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Incidence  Angle 

Figure  6.2  Comparison  of  small  perturbation  method  with  measured 
backscatter  at  Ku  band  from  CRRELEX'90.  In  this  case  the  correlation  length 
was  adjusted  to  achieve  a  better  fit.  Theory  for  the  cross-polarized  response  is 
the  integral  equation  method.  All  theoretical  curves  are  drawn  as  lines,  with 
the  measurements  as  various  symbols.  This  demonstrates  that  the  primary 
backscatter  mechanism  is  surface  scattering.  Values  used  in  model  are 
oh=1.36  mm,  Lc=3.52  cm  and  er=3.75-j0.2. 


Figure  6.3  shows  the  C-band  results,  using  the  computed  surface  parameters.  Figure 
6.4  shows  the  theoretical  results  obtained  using  the  correlation  length  that  best  fit  the 
Ku-band  results.  The  theory  and  experiment  generally  match  at  both  frequencies  for 
the  co-polarized  response.  The  measured  eft  and  that  calculated  from  the  small 
perturbation  method  agree  to  within  ±  3  dB  in  terms  of  mean-square  error.  However, 
the  theory  and  measurements  only  agree  grossly  for  the  cross-polarized  response, 
especially  at  C-band.  The  measured  return  at  low  angles  of  incidence  is 
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Incidence  Angle 

Figure  6.3  Comparison  of  small  perturbation  method  with  measured 
backscatter  at  C  band  from  CRRELEX'90.  Theory  for  the  cross-polarized 
response  is  the  integral  equation  method.  All  theoretical  curves  are  drawn  as 
lines,  with  the  measurements  as  various  symbols.  This  demonstrates  that  the 
primary  backscatter  mechanism  is  surface  scattering.  Values  used  in  model 
are  Oh=l-36  mm,  Lc=17.6  mm  and  Ef=3.75-j0.2. 


underestimated  by  the  theory  by  as  much  as  15  dB,  while  at  higher  angles  the  return  is 
overestimated  by  up  to  8  dB  at  C  band.  The  cross-polarized  measurements  and  the 
theory  only  agree  on  the  general  level  of  the  backscatter,  but  the  angular  response  of 
the  cross-polarized  measurements  is  not  predicted  by  the  method  shown  in  section  6.3 
with  the  surface  parameters  selected. 

In  both  cases  the  coherent  component  is  plotted  near  nadir  to  show  that  the  increase  in 
backscatter  near  nadir  is  due  to  the  coherent  component.  From  the  limits  imposed  on 
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Small  Perturbation 


Incidence  Angle 

Figure  6.4  Comparison  of  small  perturbation  method  with  measured  backscatter  at  C 
band  from  CRRELEX'90.  In  this  case,  the  correlation  length  was  used  from  the  better 
fit  obtained  at  Ku-band.  Theory  for  the  cross-polarized  response  is  the  integral 
equation  method.  All  theoretical  curves  are  drawn  as  lines,  with  the  measurements  as 
various  symbols.  This  demonstrates  that  the  primary  backscatter  mechanism  is 
surface  scattering.  Values  used  in  model  are  oh=1.36  mm,  Lc=35.2  mm  and  £r=3.75- 
j0.2. 

the  surface  parameters  for  the  small  perturbation  method  (Table  6.2)  we  see  that  all 
criteria  are  met  for  the  C-band  case.  In  the  Ku-band  case  the  surface  correlation 
length  is  slightly  longer  than  that  allowed  by  these  restrictions.  It  must  be 
emphasized  that  these  restrictions  are  general  guidelines  and  were  brought  about  by 
numerical  experimentation  [Chen,  1988].  In  fact  in  some  cases  the  limit  on 
correlation  length  is  considered  to  be  kg*L<6  [Ulaby,  et,  al,  1986]. 
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For  computation  of  error  measures  the  combination  of  coherent  and  incoherent 
backscatter  are  plotted  versus  the  measurements  in  figures  6.5  through  6.8.  These 
correspond  to  figures  6. 1  through  6.4,  but  have  the  result  of  adding  the  coherent  and 
incoherent  components.  The  resulting  mean-square  error  between  the  measurements 
and  the  theoretical  results  are  shown  in  table  6.4.  These  comparisons  demonstrate 
that  the  primary  backscatter  mechanism  for  this  saline  ice  is  the  scattering  from  the 
surface.  However,  the  theory  and  experiment  are  far  from  a  perfect  match.  This  may 
be  due  to  underestimating  the  surface  correlation  length  from  photographs  of  cross 
sections  of  samples  from  the  ice  sheet  or  changes  in  the  surface  roughness  across  the 
ice  sheet.  Another  cause  for  the  slight  discrepancy  between  theory  and  measurement 
at  Ku-band  is  that  these  data  are  on  the  edge  of  the  limits  for  applying  small 
perturbation  theory.  To  avoid  these  kinds  of  limitations  integral  equation  methods 
may  be  more  applicable  [Fung  and  Pan,  1987]. 


Incidence  Angle 

Figure  6.5  Sum  of  coherent  and  incoherent  theoretical  backscatter  components 
compared  to  measurements  at  Ku  band.  This  corresponds  to  the  plot  of  figure  6.1. 
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Small  Perturbation 


Incidence  Angle 

Figure  6.6  Sum  of  coherent  and  incoherent  theoretical  backscatter  components 
compared  to  measurements  at  Ku  band.  This  corresponds  to  the  plot  of  figure  6.2. 


Incidence  Angle 

Figure  6.7  Sum  of  coherent  and  incoherent  theoretical  backscatter  components 
compared  to  measurements  at  C  band.  This  corresponds  to  the  plot  of  figure  6.3. 
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Small  Perturbation 
Parametric  Fit 

Mar.  11  C  band 


Incidence  Angle 

Figure  6.8  Sum  of  coherent  and  incoherent  theoretical  backscatter  components 
compared  to  measurements  at  C  band.  This  corresponds  to  the  plot  of  figure  6.4. 


Table  6.4.  Mean-square  error  measures  (in  dB)  of  comparisons  between 

theory  and  measurement 

C-band 

ijjggpi 

Ku-band 

1 

VV 

3.8 

3.9 

2.7 

2.8 

HH 

2.7 

3.7 

2.4 

3.9 

VH 

8.5 

7.1 

9.8 

10.7 

HV 

10.9 

8.2 

8.7 

9.9 

VH  (15°-55°) 

8.4 

5.4 

7.3 

3.0 

HH  (15°-55°) 

11.9 

8.4 

6.0 

2.3 

WEMEM 

7.4, 7.8 

7.1, 7.5 

For  the  January  portion  of  the  experiment,  the  surface  parameters  were  computed 
from  a  sample  taken  on  20  January.  The  closest  time  frame  in  which  backscatter 
measurements  were  performed  was  on  24  January.  There  were  no  C-band 
measurements  made  in  late  January,  but  both  C-band  and  Ku-band  measurements 
were  made  on  9  January.  Figure  6.9  shows  the  Ku-band  Kirchoff  model  and  small 
perturbation  method  results  compared  to  backscatter  data  from  24  January.  This 
shows  that  the  Kirchoff  approach  may  be  slightly  more  appropriate  than  the  small 
perturbation  method,  but  the  characteristics  of  this  surface  lies  outside  the  regions  of 
validity  for  both  theories.  This  violation  is  not  severe,  though,  and  the  two  theories 
are  fairly  close  for  large  angles.  This  shows  that  in  this  case  also,  the  surface  scatter 
dominates  the  co-polarized  return.  Since  there  were  no  data  taken  on  20  January,  data 
from  24  January  are  shown  as  the  measured  data.  The  time  between  the  backscatter 
measurement  and  sampling  for  surface  parameterization  may  be  the  primary  reason 
for  the  difference  between  measurements  and  theory  for  this  case. 
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CRRELEX  '90 
Kirchoff  Model  Results 


Incidence  Angle 

Figure  6.9  Comparison  of  Ku-band  measurements  from  24  January,  with 
Kirchoff  model  and  small  perturbation  method.  Surface  parameters  were 
characterized  from  slab  taken  on  20  January.  These  parameters  are  ah=2.16 
mm,  Lc=  4.7  mm  and  6r=3.75-j0.2. 
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7.  CRRELEX'90  Conclusions 


7.1.  Conclusions 

Measurements  of  microwave  backscatter  from  simulated  sea  ice  demonstrate  that  the 
dominant  backscatter  mechanism  is  the  surface  scatter  when  an  exponential 
correlation  function  is  assumed  for  the  surface.  Both  the  co-polarized  and  cross- 
polarized  measurements  compare  favorably  with  the  predictions  of  surface  scattering 
models  at  two  frequencies.  The  surface  parameters  used  in  the  models  were  measured 
directly  from  photographs  of  cross  sections  of  the  saline  ice.  Therefore,  since  two 
frequencies  of  backscatter  measurements  were  predicted  accurately  by  surface  scatter 
models  the  predominate  backscatter  is  from  the  surface.  It  is  possible,  assuming 
Gaussian  correlation  for  the  surface,  to  come  to  the  conclusion  that  the  volume  plays 
a  significant  role  in  the  backscatter  from  saline  ice  [Coutu,  personal  communication, 
1992].  The  autocorrelation  of  sea  ice  surfaces  though  tend  to  be  closer  to  exponential 
than  Gaussian. 

These  measurements  were  carefully  calibrated  using  a  complex  vector  error 
correction  scheme  to  remove  polarization  coupling  and  other  systematic  effects.  The 
dynamic  range  of  polarimetric  measurements  made  by  scatterometers  are  limited  by 
coupling  and  polarization  leakage  in  the  antenna  system  and  RF  section.  These 
effects  must  be  accounted  for  and  corrected  if  one  is  to  obtain  accurate 
characterization  of  the  polarimetric  response  of  distributed  targets. 

This  demonstrates  the  need  for  careful  calibration  in  experiments  of  this  type, 
especially  when  cross-polarized  backscatter  measurements  are  performed.  The 
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systematic  effects  on  the  cross-polarized  backscatter  introduce  significant  error  and 
may  lead  to  misinterpretation  of  the  results.  This  is  a  vital  consideration  in  the  study 
of  radar  backscatter  from  sea  ice  using  fully  polarimetric  systems. 

7.2  Limitations  and  constraints 

There  are  limitations  to  the  accuracy  of  the  measurements  that  can  be  obtained  when 
applying  the  vector  error  calibration  technique.  These  limits  depend  on  how  carefully 
the  measurements  can  be  performed  in  a  practical  situation. 

The  signal-to-noise  ratio  (SNR)  of  each  calibration  target  measurement  must  be 
sufficiently  large.  This  is  a  problem  when  using  a  metal  sphere,  for  example,  because 
of  its  low  radar  cross  section.  An  SNR  of  less  than  10  dB  may  cause  enough  phase 
error  to  distort  the  calibration,  resulting  in  measurement  inaccuracies.  The  alignment 
of  the  calibration  targets,  especially  in  the  case  of  the  dihedral  comer  reflectors  is 
vital.  In  one  dimension,  the  rotation  of  the  diplane  will  affect  the  amount  of 
depolarization  attributed  to  the  system.  An  error  in  aligning  the  rotation  angle  may 
cause  a  mistake  in  characterization  of  the  system's  polarization  response.  Any  pitch 
of  the  diplane,  on  the  other  hand  will  simply  cause  a  reduction  in  the  absolute  level  of 
the  return  power.  This  effect  is  taken  into  account  in  the  normalization  of  the  power 
returns  before  computing  the  calibration  matrix. 

Since  wide-beam  antennas  were  used  for  this  experiment  their  effect  on  the 
measurement  must  be  considered.  In  using  the  wide-beam  correction  scheme 
mentioned  in  section  4.5  and  described  in  [Wang  and  Gogineni,  1990],  these  antenna 
effects  are  reduced.  The  wide-beam  antenna  actually  helps  in  practical 
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implementation  of  the  vector  calibration  because  of  the  reduced  sensitivity  with 
respect  to  alignment  than  for  narrow-beam  antenna  systems. 

The  effect  of  the  limitations  due  to  rotation  and  pitch  of  the  calibration  targets  as  well 
as  random  errors  to  low  SNR  is  to  cause  some  limit  on  the  minimum  measurable 
scattering  coefficient.  This  dynamic  range  is  significantly  larger  than  when  applying 
the  simple  scalar  calibration  as  evidenced  by  the  increase  in  dynamic  range  of  the  o° 
measurements  shown  in  sections  5  and  6. 


87 


88 


References 

Arcone,  S.A.,  A .J:  Gow,  and  S.  McGrew,  (1986),  "Microwave  Dielectric,  Structural 
and  Salinity  Properties  of  Simulated  Sea  Ice,"  IEEE  Transactions  on 
Geoscience  and  Remote  Sensing,  vol.  24,  no.  6,  pp.  832-839. 

Balanis,  C.A.,  (1982),  Antenna  Theory,  New  York  :  Harper  &  Row  Publishers,  pp. 
532-587. 

Beaven,  S.G.,  S.P.  Gogineni,  S.  Tjuatja,  A.K.  Fung,  and  J.  Bredow  (1991)  "Radar 

Backscatter  from  Saline  Ice:  Measurements  and  Comparison  of  Measurements 
with  Theory,"  IGARSS'91,  vol.  0,  pp.  425-428. 

Bredow,  J.W.,  (1989),  "A  laboratory  Investigation  Into  Microwave  Backscattering 
From  Sea  Ice,"  Ph.D.  Dissertation,  University  of  Kansas. 

Chen,  T.,  T.  Chu,  and  F.  Chen,  (1991),  “A  NeBand  Polarimetric  Measurement 

System,”  IEEE  Transactions  on  Antennas  and  Propagation,  vol.  39,  no.  8, 
pp.  1188-1192. 

Chen,  M.F.,  and  A.K.  Fung,  (1988),  "A  Numerical  Study  of  the  Regions  of  Validity 
of  the  Kirchoff  and  Small-Perturbation  Rough  Surface  Scattering  Models," 
Radio  Science,  Vol.  23,  No.  2,  pp.  163-170,  March-April,  1988. 

Coutu,  Pierre,  Massachusetts  Institute  of  Technology,  personal  communication,  Dec. 
1992. 


89 


Drinkwater,  M.R.,  R.  Kwok,  D.  Winebrenner,  and  E.  Rignot,  (1991),  "Multi- 
frequency  Polarimetric  SAR  Observations  of  Sea  Ice,"  Journal  of 
Geophysical  Research,  vol.  96  no.  Cll,  pp.  20679-20698. 

Fung,  A.K.  and  H.J.  Eom,  (1983),  "Coherent  scattering  of  of  a  Spherical  Wave  from 
an  irregular  Surface,"  IEEE  Transactions  on  Antennas  and  Propagation,  vol. 
TAP-31,  pp.  68-72. 

Fung,  A.K.,  and  G.W.  Pan,  (1987),  "A  scattering  model  for  perfectly  conducting 
random  surfaces  I.  Model  development,"  International  Journal  of  Remote 
Sensing,  Vol.  8,  no.  11,  pp.  1579-1593. 

Houghton,  J.T.,  G.J.  Jenkins,  and  J.J.  Ephraums,  (1990),  "Climate  change  -  the 
IPCC  scientific  assessment,  Cambridge  University  Press. 

Iizuka,  K.,  and  A.P.  Freundorfer,  (1984),  “Step-frequency  radar”.  Journal  of  Applied 
Physics,  vol.  56,  no.  9,  pp.  2572-2583. 

Kong,  J.A.,  (1990),  Electromagnetic  Wave  Theory.  Wiley-Interscience,  New  York. 

Kwok,  R.,  M.  Drinkwater,  A.  Pang,  and  E.  Rignot,  (1991),  "Characterization  and 
Classification  of  Sea  Ice  in  Polarimetric  SAR  Data,"  IGARSS'91,  vol.I,  pp. 
81-84. 


Love,  A.W.,  (1962),  "The  Diagonal  Horn  Antenna,"  The  Microwave  Journal ,  pp. 
117-122,  March,  1962. 


90 


Parkinson,  C.L.,  J.C.  Comiso,  HJ.  Zwally,  D.J.  Cavalieri,  P.  Gloerson,  and  W.J. 
Campbell,  (1987)  Arctic  Sea  Ice.  1973-1976:  Satellite  Passive-Microwave 
Observations.  Washington  DC:  NASA  Scientific  and  Technical  Information 
Branch. 

Rice,  S.O.  (1951),  "Reflection  of  Electromagnetic  Waves  from  Slightly  Rough 
Surfaces,  Communications  in  Pure  and  Applied  Mathematics,  vol.  4,  pp. 
361-378. 

Riegger,  S.,  and  W.  Wiesbeck,  (1989),  "Wide-band  Polarimetry  and  Complex  Radar 
Cross  Section  Signatures,"  Proceedings  of  the  IEEE,  vol.  77,  no.  5,  pp.  649- 
658. 

Swift,  C.T.,  K.  St.  Germain,  K.  Jezek,  S.P.  Gogineni,  A.  Gow,  D.  Perovich,  T. 

Grenfell,  R.  Onstott,  (1992),  Laboratory  Investigations  of  the  Electromagnetic 
Properties  of  Artificial  Sea  Ice,  in:  Microwave  Remote  Sensing  of  Sea  Ice, 
F.D.  Carsey  (ed.),  American  Geophysical  Union,  Monograph  68,  pp.  177-200. 

Ulaby,  F.T.,  R.K.  Moore,  and  A.  K.  Fung,  (1981),  Microwave  Remote  Sensing: 
Active  and  Passive,  vol.  I,  Norwood,  MA:  Artech  House. 


91 


Ulaby,  F.T.,  R.K.  Moore,  and  A.  K.  Fung,  (1982),  Microwave  Remote  Sensing: 
Active  and  Passive,  vol.  II,  Norwood,  MA:  Artech  House. 

Ulaby,  F.T.,  R.K.  Moore,  and  A.  K.  Fung,  (1986),  Microwave  Remote  Sensing: 
Active  and  Passive,  vol.  IE,  Norwood,  MA:  Artech  House. 

Ulaby,  F.T.  and  C.  Elachi,  (1990),  Radar  Polarimetrv  for  Geoscience  Applications. 
Artech  House,  Norwood,  MA. 

Valenzuela,  G.R.  (1967),  "Depolarization  of  EM  Waves  by  slightly  Rough  Surfaces," 
IEEE  Transaction  on  Antennas  and  Propagation,  vol.  TAP-15,  pp.  552-557. 

Wadhams,  P.,  (1992)  "Remote  Sensing  of  Snow  and  Ice,"  presented  at  University  of 
Dundee  Remote  Sesing  Summer  School:  Remote  Sensing  and  Global  Climate 
Change,  Unversity  of  Dundee,  Scotland,  UK,  July,  1992. 

Wang,  Q.,  and  S.P.  Gogineni,  (1991),  “A  Numerical  Procedure  for  Recovering 
Scattering  Coefficients  from  Measurements  with  Wide-Beam  Antennas”, 
IEEE  Trans,  on  Geoscience  and  Remote  Sensing,  vol.  25,  no.  5,  pp.  778-783. 

Wiesbeck.W.,  and  S.  Riegger,  (1991),  “A  Complete  Error  Model  for  Free-Space 
Polarimetric  Measurements,”  IEEE  Transactions  on  Antennas  and 
Propagation,  vol.  39,  no.  8,  pp.  1 105- 1111. 


92 


Wu,  S.T.,  and  A.K.  Fung,  (1972),  "A  Noncoherent  Model  for  Microwave  Emmisions 
and  BacksCattering  from  the  Sea  Surface,"  Journal  of  Geophysics  Research, 
vol.  77,  pp.  5917-5929. 

Acknowledgement : 

The  authors  would  like  to  thankTerry  Tucker  and  John  Govoni  at  CRREL  for  their 
assistance  in  the  experiment.  We'd  also  like  to  thank  Jim  Ronnau  and  Jeff  Young  of 
KU  for  their  participation  in  the  experiment.  This  work  was  funded  in  part  by  an 
Office  of  Naval  Research  (ONR)  doctoral  fellowship,  as  well  as  ONR  grant 
#N00014-89-J-1456. 


94 


Appendix  A 


CRRELEX'90  Scalar  calibrated  results  of  scattering  from  simulated  sea  ice 
without  beamwidth  correction  for  both  C  band  and  Ku  band 
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CRRELEX'90  VECTOR  calibrated  results  of  scattering  from  simulated  sea 
ice  without  beamwidth  correction  for  both  C  band  and  Ku  band 
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CRRELEX'90  Scalar  calibrated  results  of  scattering  from  simulated  sea  ice 
with  beam  width  correction  for  both  C  band  and  Ku  band 
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Appendix  D 


CRRELEX’90  VECTOR  calibrated  results  of  scattering  from  simulated  sea 
ice  with  beamwidth  correction  for  both  C  band  and  Ku  band 
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